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摘 X: 蒸 散 发 (ET) 是 陆地 水 循环 过 程 的 重要 组 成 部 分 ， 同 时 也 是 区 域 能 量 平衡 以 及 水 量 平衡 的 关键 环节 ， 精 
确 估算 ET， 对 于 提高 水 分 利用 效率 以 及 优化 区 域 用 水 结构 具有 重要 意义 。 本 文 利 用 黑河 重大 计划 观测 数据 ， 对 
比 了 考虑 CO. 浓度 和 不 考虑 CO. 浓度 对 玉米 冠 层 影响 的 冠 层 阻力 模型 ， 分 别 将 其 耦合 到 双 源 的 
Shuttleworth-Wallace(S-W) 模 型 中 ,并 利用 这 两 种 模型 分 时 段 对 玉米 整个 生育 期 内 半 小 时 尺度 上 的 ET 进行 模拟 ， 
利用 涡 度 相关 实测 数据 对 模型 进行 验证 ， 最 后 分 别 对 影响 玉米 冠 层 阻力 的 气象 要 素 和 影响 ET 的 阻力 参数 进行 敏 
感性 分 析 , 探寻 大 气 CO 浓度 改变 条 件 下 黑河 中 游 绿 洲 区 玉米 不 同 生长 阶段 的 农田 耗 水 规律 .结果 表明 : 本 文 所 
修正 的 考虑 CO 浓度 对 玉米 冠 层 影响 的 冠 层 阻力 模型 耦合 到 S-W 模型 后 ， 能 够 较 精 准 地 模拟 玉米 整个 生育 期 不 
同 生长 阶段 半 小 时 尺度 上 农田 耗 水 过 程 。 敏 感性 分 析 表 明 : 各 生长 阶段 冠 层 阻力 (x ) 和 冠 层 面 高 度 到 参考 面 高 
度 间 的 空气 动力 阻力 (x2 ) 对 ET 的 影响 最 为 强烈 其 他 阻力 参数 对 ET 的 影响 不 明显 , ET 的 变化 程度 随 着 x 和 x 
的 增 大 而 减 小 。 本 文 所 修正 的 考虑 CO 浓度 影响 的 分 时 段 双 源 模型 能 够 精准 地 模拟 玉米 整个 生育 期 各 生长 阶段 
的 ET， 可 为 种 植 结构 调整 和 土地 利用 方式 改变 以 及 CO; 浓度 变化 环境 下 的 农田 蒸 散 研究 提供 参考 。 
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Abstract: Evapotranspiration (ET) is composed of two separate processes — water loss to the atmosphere from soil surface by 
evaporation (E) and water loss to the atmosphere from plant canopy via transpiration (T). ET plays a key role in energy and 


water balance in agricultural system and is also a critical process in terrestrial hydrological cycle. Accurate estimation of ET is 
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significant for improving water use efficiency and optimizing regional water use, particularly in arid and semi-arid regions. 
Although ET models have been an important tool in understanding the regulation of ET in ecological, agricultural and envi- 
ronmental sciences, the accuracy of the models is limited by aerodynamic and canopy resistance. Numerous models have been 
developed to integrate aerodynamic and canopy resistances [e.g., Penman-Monteith (P-M) model] in simulating the processes 
of response of ET, but many studies have suggested that the P-M model could produce large errors under partial or sparse 
canopy conditions because it treated plant canopy and soil surface as a single entity. Next, the dual-layered Shuttle- 
worth-Wallace (S-W) model was developed to estimate ET under different conditions. In this model, the crop ET is divided 
into two components — latent heat flux from crop and that from soil. It has been tested by various surface conditions and 
widely used because of its good performance. In this study (which used maize data of three eddy covariance observations for 
the period from May through September 2012 in Heihe River Basin, an arid area in Northwestern China), two canopy resis- 
tance models coupled for maize. Two S-W models were coupled with canopy resistance models of maize taking or non-taking 
into account the effect of atmospheric CO». Then the whole maize growth period was divided into three stages, early, middle 
and late growth stages. Then maize ET on half hour scale was simulated using the two S-W models. The performances of the 
two S-W models were validated for three different growth stages using eddy covariance field-measured data. The results 
showed that simulated maize ET by the S-W model (which took into account the effect of atmospheric CO, at every growth 
stage of three different places) best agreed with field-measured eddy covariance data. Sensitivity analysis of the revised S-W 
model (taking into account the effect of atmospheric CO;) showed that maize ET was more sensitive to canopy resistance UT ) 
and aerodynamic resistance from canopy to reference surface height (rè) at different growth stages. Therefore, it is very nec- 
essary to determine resistance parameters at different growth stages taking into account the effect of atmospheric CO, when 
calculating maize ET using the revised S-W model. 


Keywords: Shuttleworth-Wallace (S-W) model; Evapotranspiration; Atmospheric CO» concentration; Resistance parameter; Maize 


植被 蒸腾 (T) 和 土壤 蒸发 (E) 所 构成 的 蒸 散发 
(ET) 是 水 循环 过 程 的 重要 组 成 部 分 , 同时 也 是 区 域 
能 量 平衡 以 及 水 量 平衡 的 关键 环节 [5 ， 精 确 估算 
蒸 散 发 ， 对 于 提高 水 分 利用 效率 以 及 优化 区 域 用 水 


物 与 外 界 的 水 气 通 量 交 换 (0 ”研究 表明 CO, 浓度 增加 
会 使 气孔 导 度 下 降 ， 作 物 蒸腾 作用 减弱 ， 而 高 浓度 
CO: 则 能 诱 使 气孔 关闭 ,进而 对 作物 蒸 散 以 及 水 分 利用 
效率 产生 较 大 影响 5”…“]。Wand SRA CO 浓度 


结构 具有 重要 意义 F”。 气孔 作为 蒸 散 发 过 程 中 植物 
与 大 气 进 行 水 碳 交 换 的 通道 ， 其 开 闭 程度 受到 植物 
生理 和 诸多 环境 要 素 的 共同 控制 ,因此 利用 模型 模 
拟 蒸 散 变化 规律 及 其 特征 成 为 蒸 散 研究 的 有 效 手 
段 。 目 前 蒸 散 模型 主要 分 为 单 源 模型 、 双 源 模 型 和 
多 源 模型 ， 这 些 模 型 分 别 引 入 各 种 阻力 参数 来 描述 
土壤 蒸发 和 植被 蒸腾 ， 其 中 描述 了 水 分 通过 作物 冠 
层 所 需要 克服 的 阻力 一 一 冠 层 阻力 ， 其 模拟 精度 决 
XE ZK BIB RE”! 

基于 能 量 平衡 的 单 源 Penman-Monteith (P-M) 模 型 
将 作物 冠 层 和 地 表 作 为 一 个 整体 , 利用 部 分 气象 要 素 
和 作物 属性 就 可 以 估算 作物 蒸 散 量 已。 众多 研究 表明 基 
于 “大 叶 ” 假 设 的 P-M 模型 对 作物 冠 层 下 表面 性 质 具 有 
较 强 的 敏感 性 “2 ， 对 作物 完全 覆盖 地 表 时 具有 较 高 的 
精度 ,而 在 作物 初始 阶段 会 产生 较 大 误差 1。 1985 
年 Shuttleworth 和 Wallace WEEE FAA BH 
行 研究 , 引入 了 冠 层 阻力 参数 和 土壤 阻力 参数 ,建立 了 
由 作物 冠 层 以 及 冠 层 下 表面 两 部 分 组 成 的 双 源 蒸 散 模 
型 ， 该 模型 由 于 考虑 了 土壤 蒸发 ， 机 理 更 加 明确 , 因而 
1S8 T T Ze, BEI NRRL HAS BA 
气 CO, 浓 度 的 影响 , 而 大 气 CO, 浓 度 变 化 直接 影响 到 作 


加 倍 时 ，C4 作物 的 气孔 导 度 降低 29%; Morison S&P 
CO, 浓 度 升 高 情境 下 植物 的 响应 规律 进行 了 研究 , 结果 
表明 CO, 浓 度 加 信 会 使 大 田 作 物 的 气孔 导 度 降低 40%。 
在 以 气温 上 升 和 大 气 CO. 浓度 升 高 为 主要 特征 
的 全 球 变化 背景 下 ,大 气 CO, 浓度 升 高 和 水 资源 短 
缺 已 经 对 区 域 农 业 的 可 持续 发 展 产生 了 严重 影响 7。 
而 目前 气候 变化 环境 下 的 农田 蒸 散 研究 并 未 考虑 到 
KA CO, 浓度 与 其 他 环境 因子 所 产生 的 协同 效应 ， 
因此 考虑 大 气 CO. 浓度 对 作物 蒸 散 的 影响 ， 对 于 揭 
示 作 物耗 水 规律 和 提高 水 分 利用 效率 具有 重要 意义 ， 
尤其 是 在 水 分 限制 区 域 。 现 有 蒸 散 模型 中 的 冠 层 阻 
力 参数 并 未 考虑 大 气 CO. 浓度 的 影响 , 但 从 模型 长 
期 预测 的 角度 考虑 ,尤其 是 在 未 来 气候 和 生态 系统 
水 碳 平 衡 变化 的 模拟 预测 中 中， 考虑 大 气 CO, 浓度 
变化 对 冠 层 阻力 和 蒸 散 的 影响 显得 尤为 重要 。 此 外 ， 
作物 在 不 同 生育 阶段 具有 不 同 的 蒸 散 规 律 ,利用 模 
型 估算 整个 生育 期 的 蒸 散 量 ， 使 用 一 套 相 应 的 参数 ， 
不 尽 合理 ， 因 此 有 必要 结合 作物 不 同 生育 期 耗 水 规 
律 ， 分 阶段 估算 作物 蒸 散 量 ， 进 而 提高 模拟 精度 。 
本 研究 基于 Morison 等 C0 和 Easterling SO" X F 
CO; 浓度 变化 对 气孔 导 度 的 影响 规律 研究 ,利用 “ 黑 
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河流 域 生 态 -水 文 过 程 集成 研究 ”重大 计划 已 有 观测 
数据 ,以 Shuttleworth 和 Wallace 所 构建 的 双 源 模型 
(S-W) 为 基础 ， 将 玉米 (Zea mays) 生 育 期 分 为 3 SE 
要 生长 阶段 ， 对 比分 析 了 考虑 和 未 考虑 CO, 浓度 对 
黑河 中 游 绿 洲 区 玉米 冠 层 阻力 的 影响 ,并 将 其 耦合 
到 双 源 的 S-W 模型 中 ， 分 阶段 估算 玉米 蒸腾 和 土壤 
蒸发 ， 改 进 部 分 阻力 参数 ， 利 用 涡 度 相关 实测 数据 
对 模型 进行 验证 ， 筛 选 出 能 够 反映 大 气 CO, KES 
化 情境 下 的 冠 层 阻力 模型 ， 并 将 其 耦合 到 S-W 模型 
中 ,探寻 大 气 CO. 浓度 改变 条 件 下 黑河 中 游 绿 洲 区 
玉米 不 同 生长 阶段 的 农田 耗 水 规律 。 
1 研究 区 概况 与 研究 方法 
1.1 ”研究 区 概况 

黑河 中 游 绿洲 区 位 于 河西 走廊 中 部 (图 D), 介 于 
98°57'~100°S1'E, 38°32'~39°53'N, BRA 
1 131-2 891 m, RTH USE Zi HL ETT CA 
下 至 正义 峡 之 间 的 地 势 平坦 区 域 。 本 研究 区 为 典型 
的 温带 大 陆 性 气候 ， 流 域 气候 干燥 ， 降 水 稀少 ,年 
降水 量 116.8 mm, 年 蒸发 量 2 365.6 mm, KAKA 
源 分 配 严重 不 均 。 地 带 性 土壤 为 厌 棕 荒漠 土 和 灰 漠 土 ， 
非 地 带 性 土壤 主要 包括 盐 土 、 潮 土 、 风 沙土 等 。 玉 米 
是 该 地 区 最 主要 的 农作物 之 一 ， 其 生长 主要 依赖 于 


[X 3 Maize I 


地 下 水 以 及 黑河 水 灌溉 。 
12 ”数据 来 源 及 处 理 

本 研究 所 使 用 的 实测 数据 包括 : 自动 气象 站 资 
料 、 通 量 数据 、 作 物 株 高 等 ， 该 数据 均 来 自 黑河 计 
划 数 据 管 理 中 心 (http://heihedata.org/)。 自 动 气象 站 
观测 气象 要 素 包 括 风 速 、 风 向 、 气 温 、 降 水 量 、 湿 
度 等 。 叶 面积 指数 (LAT) 来 源 于 美国 国家 航天 局 
(NASA) 发 布 的 MCDISA3H 叶 面 积 指 数 (LAD) 产 品 
(https://search.earthdata.nasa.gov/)， 该 数据 为 4 d 合 
成 的 500 m 分辩 率 L4 级 叶 面 积 指 数 产 品 ,， 采 用 3 次 
样 条 插值 将 其 均匀 内 插 到 日 尺度 ， 并 假设 一 天 内 玉 
米 的 叶 面 积 指 数 不 发 生变 化 。 本 研究 使 用 通 量 观测 算 
阵 中 Site 8(100°22'35"E, 38?5221"N; 1 550.06 m) 的 数 
据 率 定 模型 参数 ，Site 11(100?20/31"E, 38?52'12"N; 
1575.65 m). X35835(100922'20"E, 38951'20"N; 
1 556.06 m) 的 数据 用 于 模型 的 验证 (图 1)。3 个 站 点 观 
测 时 间 序 列 为 整个 玉米 生长 季 (2012 年 5 月 初 到 9 月 
下 旬 )， 观 测 时 间 为 9:00 一 18:30。3 站 点 的 通 量 观测 系 
统 主要 由 三 维 超声 风 温 仪 (CSAT3, Campbell Scientific, 
USA) 和 开路 COJH;,O 红外 气体 分 析 仪 Li-7500A， 
Li-Cor Imc.，USA) 组 成 ， 原 始 数据 采样 频率 为 10 Hz, 
每 30 min 输出 1 组 平均 通 量 值 。 通 量 站 点 下 垫 面 作 


站 点 Experimental site 


| 积 雪 Snow cover — 黑河 Heihe River 


Eg HAE Alfalfa 
[Rest Grass land BR 其 他 农 地 Other farm land 
D iis Broad leaved forest J) 3 Barley 

Ett Bare land 
[X 4n Coniferous forest B 湿地 Wet land 


= 建设 用 地 Construction land 


B 水 体 Water body 


棉花 Cotton 


[9] iE HG Bush fallow 


1 黑河 中 游 绿洲 区 土地 利用 类 型 及 试验 站 点 分 布 


Fig.1 Land use types of Heihe River basin and distribution of the experiment sites 
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物 均 为 玉米 ， 垄 距 为 50.8 cm, 行距 为 43.3 cm， 株 距 
为 22 cm, 玉米 株 高 由 实测 数据 采用 3 次 样 条 插值 均 
匀 内 插 到 半 小 时 尺度 。 原 始 涡 度 相 关 仪 器 信息 及 数据 
处 理 参考 文献 [30]， 原始 数据 经 过 了 野 点 值 噜 除 、 延 迟 
时 间 校 正 、 坐 标 旋转 、 角 度 订 正 以 及 严格 的 质量 控制 
等 步骤 , 但 还 需要 对 部 分 数据 进行 再 处 理 : 1) 能 量 闭合 
度 检查 , 将 能 量 闭 合 度 超出 0.5~1.5 的 数据 剔除 ， 能量 
未 闭合 的 数据 采用 文献 [31] 中 的 方法 进行 强制 闭合 ; 
2) 别 除 空 值 数 据 和 异常 数据 ， 如 实测 的 潜 热 通 量 
(LE)«0 W:m^, r5,,»2 000 sm 或 euo<0 sm! (r^ A 
利用 Penman-Monteith 公式 反 推 的 冠 层 阻力 值 ， 本 
文 定 义 为 实测 值 )。 

本 研究 将 玉米 整个 生育 期 分 为 3 个 生长 阶段 : 
生育 前 期 (2012 年 6 月 6 日 一 7 月 20 日 , 播种 期 一 出 
苗 期 一 五 叶 期 一 拔节 期 一 抽穗 期 一 吐 丝 期 )， 该 阶段 
玉米 快速 生长 , LALA AK, HRE RANDMA 
到 稠密 ， 从 土壤 蒸发 大 于 蒸腾 向 蒸腾 逐渐 占 主 导 地 
位 过 渡 ; 生育 中 期 (2012 年 7 月 21 日 一 8 月 31 日 , 吐 
丝 期 一 成 熟 期 )， 该 阶段 地 表 逐 渐 全 部 覆盖 ,玉米 蒸 
腾 逐 渐 占 主导 地 位 ,土壤 蒸 发 逐渐 减弱 ; 生育 后 期 
(2012 年 9 月 1 一 20 日 , 成 熟 期 一 收获 期 )， 此 阶段 玉 
米 基本 成 熟 ， 叶片 逐 渐 豪 老 ， 叶片 气孔 导 度 减弱 ， 
群体 冠 层 阻力 增 大 ， 作 物 蒸腾 明显 减弱 ,土壤 蒸发 
较 上 一 阶段 有 所 增强 。 

13 ”模型 及 参数 确定 
153.4. S-W 模型 

Shuttleworth 和 Wallace!f£ 1985 年 以 P-M 模型 
为 基础 ， 对 稀疏 覆盖 表面 下 的 菊 散 进行 研究 ， 引 入 
冠 层 阻力 和 土壤 阻力 参数 ,建立 了 由 作物 冠 层 以 及 
冠 层 下 部 所 组 成 的 双 源 蒸 散 模 型 。 模 型 形式 如 下 : 

和 ET 2 AT +AE = C,PM, + C,PM, (1) 


PM. = M «| (pC, VPD - Ari 4, )/ (o? a à 
petite 


Aye iz +r ] 


1 


*o 1+[RR,/R, (R. +R, )] 7 
1 

”1+[RR /R,(R, +R, )] > 

R, =(A+y)r° tyre (6) 

R, -(^*«y)r «r5 (?) 

R, -(A*y)n (8) 


林 等 : 分 时 段 修正 双 源 模型 在 西北 干旱 区 玉米 蒸 散 量 模拟 中 的 应 用 Fr hinay jy 81637 
A=R,-G (9) 
A, = Ras -G (10) 
Ris x R,exp(-cLAI) (11) 


式 中 : 和 ET 为 总 的 潜 热 通 量 , 包括 作 物 蒸腾 (和 7 ) 和 
土壤 蒸发 (XE ); PM.、 PM 分别 为 冠 层 潜 热 通 量 和 
冠 层 下 部 潜 热 通 量 ，W:m ”， 分 别 用 于 计算 作物 蒸腾 
MERRE; Co C, 分 别 为 冠 层 阻力 分 配 系数 和 土 
塘 表 面 阻力 分 配 系数 ; 心 为 冠 层 面 高 度 到 参考 面 高 
度 间 的 空气 动力 阻力 ， sm -; ww 为 地 表面 到 冠 层面 
839272777 [H72, sm 5; ww 为 冠 层 叶 面 边界 层 空 气 
动力 阻力 , sm”; 7? 为 地 表面 土壤 阻力 , sm’; x 为 冠 层 
阻力 , sm '; A 为 饱和 水 汽 压 随 温度 的 变化 率 , kPa C; 
7 为 干 湿 表 常数 , kPa C; o 为 空气 密度 , kem ^; CG, 为 
空气 定 压 比 热 , JkPa C; VPD 为 饱和 水 汽 压 差 , kPa; 
A 为 总 有 效能 量 , Wm”; 4, 为 地 面 有 效能 量 , Wm”; 
Ry Ros 分 别 为 冠 层 上 方 和 地 表面 所 接收 的 太阳 净 辐 
Bj, Wm”; G 为 土壤 热 通 量 , W-m 7; c 为 消光 系数 ， 本 
研究 参考 文献 [4,8] 采 用 线性 插值 获得 ,LAI 为 叶 面积 
指数 , mm” 
1.3.2 ”各 阻力 参数 的 确定 

冠 层 面 高 度 到 参考 面 高 度 间 的 空气 动力 阻力 x* 
和 地 表面 到 冠 层面 的 空气 动力 阻力 x 分别 由 下 式 计 
BRB oe, 


r= : o 2) k 二 全 天 上 (12) 
kus h, —d nK, h, 
s hexp(n) —nZo zotd 
= = 13 
ra am ed e 


式 中 :大 为 卡 曼 常 数 , 0.41; u, 为 摩擦 速度 , ms; d AS 
平面 位 移 , m; h 为 冠 层 平均 高 度 , m; z 为 参考 高 度 , m; 
Z 为 粗糙 长 度 , m; 2Z4 =0.01 m; Ky AEA mI BA 
ZR, ms; n ABARA WAI BOISE, HX 
2.6P77- 其 他 参数 计算 参考 文献 [33]。 

冠 层 叶 面 边界 层 空 气动 力 阻力 x* 和 地 表面 土壤 
阻力 必 E FERRI 


Ty 


12 14 
^ 2LAI (14) 
HS = Kind (0) (15) 
0 
reas ss (16) 


式 中 : n 为 边界 层 平均 阻力 ， sm 本文 取 50 sm |), 
LAI AHMAR, mm’; 5. 29-3 E E ER 


> smin 


力 , «sm !, ASH 100 sm 310 为 0~100 cm 深 
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度 内 平均 土壤 含水 量 , cm*:cm 777, 0,73 0-100 cm 
深度 内 平均 田间 持 水 量 , cm cm >。 

考虑 叶 面 积 指 数 LAI、 净 辐射 R HRERS 
水 量 96、 饱 和 水 汽 压 着 VPD、 气 温 7, 以 及 大 气 CO, 
浓度 等 要 素 对 冠 层 阻 力 的 影响 ,构建 了 考虑 CO, 
度 影 响 (x1 ) 和 未 考虑 CO. 浓度 (7 ) 影 响 的 冠 层 阻力 
模型 ， 具 体 如 下 B8391: 
ri = Kein {LAL f (R, )x f(6)x f(VPD)x f (T, )x f(CO,)} 


(17) 
15 = Kain /{LAIX f (R, )x f(6)x f(VPD) x f(7,)} (18) 
R,(1100+a) 
= n 1 
f(&) 1100(R, +a) d 
8-8 
8)- * 20 
f(9) 8. à, (20) 
f (VPD) = exp(-b x VPD) (21) 
S(T) =T, ^c (22) 
CO, 
f(CO,)=1+(1- 535) AgiCO, (23) 


式 中 : on GS RMF CO 浓度 和 未 考虑 CO 浓度 
影响 的 冠 层 阻力 ;AR )、 f(0). f(VPD). FT) 
f(CO,) 分 别 为 净 辐 射 、 土 壤 有 效 含 水 量 、 饱 和 水 汽 压 
Æ, 气温 以 及 CO; REWER, 本 研究 在 玉米 生育 前 
HAFIR RREN -e 为 20 sm!) 后 期 取 294 sm! CO, 
浓度 胁迫 画 数 中 Ag,CO, 为 CO; 浓度 增加 1 倍 时 ， 叶 片 
气孔 导 度 减 小 的 倍数 P， 本 研究 取 0.3; 田间 持 水 量 0. 
BY 0.34 cm^cm^, BAW 0, 取 0.1 cm^cm 5; a. by c 
均 为 待 率 定 的 经 验 系 数 。 
14 ”模型 评价 

本 文采 用 决定 系数 (R*"), 均值 偏 移 误 差 (mean 
bias error，MBE) 以 及 均 方 根 误差 (root mean square 
error，RMSE) 来 评价 模型 模拟 值 (z,) 与 实测 值 (O,,) 
之 间 的 差异 , 检验 模型 模拟 精度 。 


5 " (O, i -E m)? 
R =1 Dani Oa (24) 


Y, [oni Om + Em - Ou | 


MBE = CA - Os) (25) 
na 


RMSE - I S5, =O.) (26) 
n iz 


式 中 : 6, 为 模拟 值 , Oo 为 实测 值 ，O,, 为 实测 值 的 平 
均值 , ”为 样本 容量 。 本 研究 中 实测 值 (ET。u) 为 涡 度 
相关 所 观测 的 潜 热 通 量 值 ， 模 拟 值 为 利用 模型 模拟 
的 ET 值 (ET。 


2 模型 验证 与 分 析 
2.1 ”模型 参数 率 定 及 验证 

参数 率 定 采用 Site 8 的 通 量 数据 (2012 年 6 月 6 
日 一 9 月 208), 利用 最 小 二 乘法 对 生育 前 期 (6 月 6 
日 一 7 月 20 日 )、 生 育 中 期 (7 月 21 日 一 8 月 31 A) 
生育 后 期 (9 月 1—20 日 ) 进 行 非 线性 回归 ， 率 定 模型 
参数 ， 模 型 参数 优选 值 见 表 1, 模型 对 上 比 结果 见 图 2 
和 图 3, 采用 Site 11 和 大 满 站 的 数据 检验 模型 模拟 
精度 。 


表 1 考虑 大 气 CO,s 浓 度 (点 ) 和 未 考虑 CO, 浓 度 ( 蕊 ) 影 响 
的 玉米 不 同 生长 阶段 冠 层 阻 力 模型 参数 最 优 值 


Table 1 Optimum parameters of the two canopy resistance 
models considering ( 75 ) and non-considering (75 ) CO» con- 


centration in three growing stages of maize 


生长 阶段 模型 b : 
Growth stage Model 

生育 前 期 r$ 2 243.410 0.480 -0.350 
Early growth stage $ 2 248.128 0.480 | —0.335 
生育 中 期 rs 7 293.761 0.551 -0.264 
Middle growth stage $ 8 193.055 0.546 -0.251 
生育 后 期 rs 407.517 0.388  -0.002 
Late growth stage c 297.432 0.343  —0.038 


由 图 2 和 图 3 可 知 ， 在 参数 率 定 期 ， 将 考虑 了 
CO; 浓度 和 未 考虑 CO. 浓度 影响 的 冠 层 阻力 模型 看 
合 到 双 源 的 S-W 模型 中 ,二 者 均 能 够 较为 准确 地 反 
映 玉米 在 不 同 生长 阶段 的 蒸 散 变 化 规律 。 但 考虑 了 
CO, 浓度 对 玉米 冠 层 阻力 影响 的 S-W 模型 (图 2e) 模 
拟 精度 更 高 ,， 各 生长 阶段 蒸 散 变化 规律 均 得 到 了 更 
为 真实 的 反映 ; 而 未 考虑 CO; 浓度 对 玉米 冠 层 阻力 
影响 的 S-W 模型 (图 3e),， 所 模拟 的 生育 中 期 蒸 散 量 
偏 高 ， 进 而 降低 了 整个 生育 期 的 模拟 效果 。 

在 生育 前 期 (图 2a 和 图 3a)， 玉 米 快速 生长 , LAI 
逐渐 增 大 ， 其 呼吸 作用 较为 强烈 ， 地 表 覆 盖 度 从 稀 
玖 到 稠密 ， 该 阶段 由 土壤 蕉 发 大 于 蒸腾 向 蔡 腾 逐渐 
占 主导 地 位 过 渡 ， 此 时 农田 蒸 散 最 为 复杂 ， 这 也 降 
低 了 许多 模型 在 此 阶段 的 模拟 效果 , 但 考虑 了 CO; 
浓度 对 玉米 冠 层 阻力 影响 的 S-W 模型 (图 2a) 能 够 从 
最 大 程度 上 反映 该 阶段 玉米 蒸 散 量 对 环境 变量 的 响 
应 过 程 。 在 玉米 生育 中 期 (图 2b)， 光 照 、 热 量 和 水 分 
充足 ， 玉 米 光 合作 用 最 为 强烈 ， 对 CO. 需求 较 大 ， 
此 阶段 地 表 逐 渐 全 部 覆盖 ， 冠 层 获 得 的 能 量 主要 用 
于 干 物质 积累 ， 该 阶段 玉米 蒸腾 逐渐 占 主 导 地 位 ， 
土壤 蒸发 进一步 减弱 ， 相 对 于 前 一 阶段 蒸 散 较 为 稳 
定 ， 对 比 图 2b 和 图 3b 可 以 看 出 ,考虑 了 CO, 浓度 对 
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图 2 考虑 CO; 浓度 影响 的 冠 层 阻力 模型 在 参 


数 率 定 期 对 不 同 生长 阶段 玉米 燕 散 量 的 模拟 效果 (Site 8) 


Fig. 2 Performance of the canopy resistance model considering CO» concentration to estimate evapotranspiration during different 
growing stages in the phase of calibration (Site 8) 
ET, 蒸 散 量 观 测 值 ; ETa: AMBER (A. ET, observed evapotranspiration; ETa: model simulated evapotranspiration. 


玉米 冠 层 影响 的 S-W 模型 能 够 很 好 地 反映 这 一 过 程 ， 
而 未 考虑 CO 浓度 对 玉米 冠 层 影响 的 S-W 模型 模拟 
效果 较 低 。 在 生育 后 期 (图 2c)， 玉 米 基 本 成 熟 ， 叶 片 
快速 衰老 叶片 气孔 导 度 减弱 ,群体 冠 层 阻力 增 大 ， 
加 之 母体 去 除 过 程 对 剩余 植株 叶片 的 损害 ， 作 物 蒸 
腾 明 显 减 弱 ， 土 壤 蒸发 较 上 一 阶段 有 所 增强 ， 相 比 
之 下 考虑 了 CO, 浓度 对 玉米 冠 层 阻 力 影响 的 S-W 模 
型 (图 2c) 仍 具有 较 高 的 模拟 精度 。 

从 图 2d. Ed 2e 和 图 3d 图 3e 可 以 看 出 ， 对 于 
Site 8 整个 生育 期 内 玉米 蒸 散 量 的 模拟 ， 考 虑 了 CO, 
浓度 对 玉米 冠 层 阻力 影响 的 S-W 模型 (图 2d) 比 未 考 
FE CO. 浓度 对 玉米 冠 层 阻力 影响 的 S-W 模型 (图 3d) 
模拟 效果 更 好 ， 二 者 模拟 值 与 实测 值 之 间 的 R 分 别 
为 0.96、0.93, 但 是 前 者 的 RMSE 和 MBE 要 比 后 者 


小 。 对 整个 生育 期 而 言 ， 前 者 仪 在 生育 后 期 (图 2e) 存 
在 小 部 分 高 估 现 象 ， 而 后 者 在 整个 生育 期 (图 3e) 均 存 


在 高 估 现 象 , 因此 将 考虑 了 CO. 浓度 影响 的 冠 层 阻 


力 模型 耦合 到 sw 模型 中 后 , 能够 更 为 准确 地 用 于 
估算 玉米 整个 生育 期 内 不 同 生 长 阶段 的 蒸 散 量 。 
2.2 ”模型 进一步 验证 

将 考虑 了 CO. 浓度 影响 和 未 考虑 CO, 浓度 影响 
的 冠 层 阻 力 模型 耦合 到 双 源 的 S- w 模型 中 ， 分 别 对 
Site 11 和 大 满 站 玉米 的 蒸 散 量 进行 模拟 ， 模 拟 效果 
见 图 4- 图 7。 由 图 4、 图 5 可 知 : Hi Site 11, 考虑 
了 CO 浓度 对 玉米 冠 层 阻力 影响 的 S-W 模型 (图 4) 


和 未 考虑 CO, 浓 度 影 响 的 S-W 模型 (图 5), 均 能 够 较 
好 地 模拟 整个 生育 期 内 不 同 生 长 阶段 玉米 的 蒸 散 


量 。 相 比 之 下 , 考虑 了 CO, 浓度 对 玉米 冠 层 阻力 影 
响 的 模型 ， 能 够 更 好 地 模拟 玉米 3 个 生长 阶段 蒸 散 
量变 化 , 3 个 阶段 模拟 值 与 实测 值 之 间 的 决定 系数 分 
ASA 0.96、0.98 和 0.88， 整 个 生育 期 (图 4d) 的 决定 
系数 达 0.96, 各 生长 阶段 模拟 值 与 实测 值 基本 吻合 ， 
误差 很 小 。 而 未 考虑 CO 浓度 对 玉米 冠 层 阻力 影响 
BY S-W 模型 ， 对 于 生育 前 期 (图 5a) 和 后 期 (图 5c) 模 


http://www.ecoagri.ac.cn 


640 


中 国生 态 农业 学 报 2017 


55-25 75 


800 


RMSE-33.76 í 
MBE=12.82 €. * 
600} R°=0.95 E 
n=620 : 
400 - 
200} m 
e. a 生育 前 期 
€ " Early growth stage 
E 0 200 400 600 
^ 600 
i RMSE=58.08 
MBE--13.63 : 
R-0.75 > 
400r 2278 
200 


c 生育 后 期 
Late growth stage 


400 


1 000 Fe MSE-130.74 


MBE-108.49 
R=0.84 
n=557 


800 


RMSE=87.60 
MBE-44.39 
R’=0.93 


800 1 000 


ol 
2012-06-05 


2012-06-22 2012-07-10 


2012-07-27 


2012-08-14 2012-08-31 2012-09-18 


日 期 (年 -月 -日 ) Date (year-month-day) 


3 ”未 考虑 CO, 浓度 影响 的 冠 层 阻力 模型 在 参 


数 率 定 期 对 不 同 生 长 阶段 玉米 燕 散 量 的 模拟 效果 (Site 8) 


Fig.3 Performance of the canopy resistance model non-considering CO; concentration to estimate evapotranspiration during dif- 
ferent growing stages in the phase of calibration (Site 8) 
ET. 蒸 散 量 观 测 值 ; ETa: AAS ERA. ET: observed evapotranspiration; ET: model simulated evapotranspiration. 


拟 效 果 较 差 ， 误 差 较 大 , 3 个 阶段 模拟 值 与 实测 值 之 
间 的 决定 系数 分 别 为 0.95、0.97 和 0.64。 对 于 整个 
育 期 (图 Se) 玉 米 蒸 散 量 的 模拟 值 均 存在 低估 现象 ， 
尤其 是 在 中 后 期 ， 这 主要 是 因为 该 阶段 光照 和 气温 
降低 ， 不 同时 刻 的 作物 蒸腾 和 土壤 蒸发 呈现 出 较 大 
的 变异 性 ， 而 未 考虑 CO, 浓度 对 玉米 冠 层 阻力 影响 
的 S-W 模型 对 这 种 变化 不 能 快速 响应 ， 使 得 部 分 模 
拟 值 偏 低 。 

同样 利用 两 种 模型 对 大 满 站 玉米 蒸 散 量 进 行 模 
拟 , 模拟 效果 见 图 6 和 图 7。 对 比 图 6 和 图 7 可 知 ， 考 
ET CO, 浓度 对 玉米 冠 层 阻力 影响 的 S-W 模型 (图 
6)， 能 够 更 好 地 模拟 大 满 站 玉米 3 个 生长 阶段 的 蒸 
散 量 , 模拟 值 与 实测 值 具有 较 高 的 一 致 性 , 3 个 阶段 
模拟 值 与 实测 值 之 间 的 决定 系数 分 别 达 0.98、0.98 
和 0.94， 整 个 生育 期 (图 6d) 的 决定 系数 达 0.97, SH 
段 的 RMSE 和 MBE 都 很 小 。 而 未 考虑 CO. 浓度 对 


玉米 冠 层 阻力 影响 的 S-W 模型 (图 7) 中 ， 对 大 满 站 玉 
米 生 育 前 期 (图 7a) 存 在 低估 现象 ， 且 不 能 很 好 地 反 
映 各 生长 阶段 玉米 蒸 散 量 的 峰值 。 由 于 后 期 土壤 蒸 
发 和 玉米 蒸腾 具有 较 大 的 变异 性 ， 使 得 后 期 的 模拟 
误差 较 大 , 3 个 阶段 模拟 值 与 实测 值 之 间 的 决定 系数 
分 别 为 0.97.0.97 和 0.80， 整 个 生育 期 (图 7d) 的 决定 
系数 达 0.95， 模 拟 精度 较 低 。 

综合 图 4. E 5. El 6 和 图 7 可 知 ,将 考虑 了 co, 
浓度 对 玉米 冠 层 阻力 影响 的 模型 耦合 到 双 源 的 S-W 
模型 后 ， 能 够 更 好 地 模拟 玉米 在 生理 结构 、 和 气象 条 
件 发 生变 化 情况 下 蒸 散 量 的 变化 ， 能 够 较为 准确 地 
区 分 土壤 蒸发 和 作物 蒸腾 ， 模 拟 值 与 实测 值 具有 较 
高 的 一 致 性 ,尤其 是 在 玉米 生育 前 期 和 中 期 ， 并 且 
模型 参数 容易 获取 且 具 有 较 好 的 适应 性 ， 这 表明 该 
模型 能 够 用 于 估算 干旱 区 覆盖 度 较 低 情 况 下 玉米 的 
AMS, DLR AAA P-M 模型 在 模拟 稀疏 覆盖 条 
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图 4 考虑 CO; 浓度 


度 影 响 的 冠 层 阻力 模型 在 模型 验证 期 对 不 同 生 长 阶段 玉米 蒸 散 量 的 模拟 效果 (Site 11) 


Fig.4 Performance of the canopy resistance model considering CO» concentration to estimate evapotranspiration during different growing 
stages in the phase of validation (Site 11) 


ET, ARB WM; ET 


件 下 的 蒸 散 所 产生 的 较 大 误差 问题 。 
2.3 ”敏感 性 分 析 

由 于 作物 自身 因素 以 及 外 部 要 素 的 改变 , 模型 
对 外 部 要 素 的 响应 存在 较 大 差异 ， 因 此 首先 分 析 
CO, 浓度 等 要 素 变 化 对 冠 层 阻力 模型 的 敏感 性 ， 其 
次 分 析 各 阻力 参数 变化 对 耦合 后 双 源 S-W 模型 的 敏 
感性 。 引 入 模型 结果 关于 变量 的 敏感 性 系数 ， 如 下 


式 所 示 : 
s, = tm (EE) Ea (27) 
“Ay, 0\ Ax, / x; ox, R 

式 中 : S, 为 模型 结果 关于 变量 x 的 敏感 系数 ， 无 量 


纲 ， 可 以 进行 模型 不 同 变量 之 间 的 敏感 性 对 比 ， 当 
S.>0 时 ,表示 模型 结果 随 变量 y 的 增 大 而 增 大 ， 
4s, <0 时 ， 表 示 模 型 结果 随 变 量 x 的 增 大 而 减 小 ， 
小 反映 了 模型 结果 随 变 量 x, 变化 的 敏感 程 
E, HARA, SE r 对 模型 结果 的 影响 越 大 。 


FR SREB. ET»: observed evapotranspiration; ET: model simulated evapotranspiration. 


首先 计算 不 同 生 长 阶段 各 要 素 变 化 +10% 和 
++30% 时 对 冠 层 阻力 的 影响 。 由 图 8 可 知 ， 在 玉米 
生育 前 期 (图 gj 和 中 期 图 8b) 对 影响 最 大 的 气象 
要 素 是 净 辐 射 R,、 土 壤 有 效 含水 量 9 和 叶 面 积 指 数 
LAI， 其 次 是 饱和 水 汽 压 差 VPD， 最 后 是 气温 Ta 和 
CO, 浓度 ， 并 且 随 着 要 素 变 幅 的 增加 而 增 大 ; 在 生 
育 后 期 (图 8c)， 玉 米 冠 层 阻 力 对 土壤 有 效 含 水量 0 
和 叶 面 积 指 数 LAI 最 为 敏感 ， 其 次 是 净 辐 射 R,。 在 
生育 中 期 (图 8b) 和 后 期 (图 8c), 4 CO. 浓度 分 别 增 
加 10% 和 30% 时 对 x 的 影响 程度 已 经 超过 
对 下 的 影响 。 

其 次 计算 玉米 不 同 生长 阶段 各 阻力 参数 变化 +10% 


了 气温 T, 


和 +30% 时 对 蒸 散 发 ET 的 影响 , 结果 见 图 9。 由 图 9 
可 以 看 出 , 在 整个 生育 期 对 ET 影响 最 大 的 阻力 参数 


是 冠 层 阻力 x*， 其 次 是 冠 层面 高 度 到 参考 面 高 度 间 的 
空气 动力 阻力 x*,，ET 对 其 他 阻力 参数 敏感 性 程度 
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5 未 考虑 CO; 浓度 影响 的 冠 层 阻 力 模型 在 模型 验证 期 对 不 同 生长 阶段 玉米 蒸 散 量 的 模拟 效果 (Site 11) 


Fig.5 Performance of the canopy resistance model non-considering CO; concentration to estimate evapotranspiration during dif- 
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6 考虑 CO; 浓度 影响 的 冠 层 阻 力 模型 在 模型 验证 期 对 不 同 生 长 阶段 玉米 燕 散 量 的 模拟 效果 (大 满 ) 
Fig. 6 Performance of the canopy resistance model considering CO; concentration to estimate evapotranspiration during different growing 
stages in the phase of validation (Site Daman) 
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Fig. 7 Performance of the canopy resistance model non-considering CO» concentration to estimate evapotranspiration during dif- 
ferent growing stages in the phase of validation (Site Daman) 
ET, KAS WIE: ET: AMS RHA. ET: observed evapotranspiration; ET,,: model simulated evapotranspiration. 


很 小 ,这 种 影响 在 生育 前 期 (图 9a) 和 后 期 (图 9b) 最 
为 强烈 。ET 的 变化 量 随 着 x* 和 的 增 大 而 减 小 , 但 
这 种 影响 作用 随 着 x* 和 x 变动 幅度 的 增 大 而 逐步 降 
低 。 分 析 可 知 ，x* 代表 了 通过 玉米 气孔 和 所 有 了 叶 表 
面 的 水 气流 所 需要 克服 的 阻力 ，x* 的 增加 意味 着 气 
流通 过 玉米 叶片 所 需要 克服 的 阻力 增 大 ，x* 的 增加 


使 得 玉米 冠 层 导 度 减 小 ， 从 而 抑制 了 水 分 流失 ,反之 
术 的 减 小 增 大 了 冠 层 导 度 ,使 得 玉米 敬 腾 更 为 旺盛 ， 
从 而 带 走 了 更 多 的 水 分 。 因 此 , 将 考虑 了 CO, 浓度 影 
响 的 冠 层 阻力 模型 耦合 到 双 源 S-W 模型 中 分 阶段 计算 
农田 蒸 散 量 时 ， 需 要 特别 注意 不 同 生 育 阶 段 对 蒸 散 量 
计算 结果 影响 较 大 的 阻力 参数 必 和 忆 的 合理 确定 。 
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图 8 玉米 生育 前 期 (a)、 中 期 (b) 和 后 期 (0) 气 象 要 素 变化 时 对 冠 层 阻力 (x ) 的 敏感 性 分 析 


Fig.8 Sensitivity analysis between climatic factors and canopy resistance ( A ) at early (a), middle (b) and late (c) 


growing stages of maize 
Ry: 冠 层 上 方太 阳 净 辐射 ; 0: 0~100 cm 深度 内 平均 土壤 含水 量 ; VPD: 饱和 水 汽 压 差 ; Ty: 气温 ; CO: AA CoO RE; LAL 叶 面 积 指数 。 Ri: net solar 
radiation above canopy; 0: average soil water content of 0—100 cm layer; VPD: vapor pressure deficit; Tą: air temperature; CO»: air CO» concentration; LAT: 
leaf area index. 


20 MERERI 40 
azar Wr Br E Sr; b F 


Um, rd us 


10 30 +10 230 2l 10 30 HO — 2430 
变化 率 Change rate (96) 
图 9 玉米 生育 前 期 (a)、 中 期 (pD) 和 后 期 (各 阻力 参数 变化 时 对 蒸 散 (ET) 的 敏感 性 分 析 


Fig.9 Sensitivity analysis between resistance parameters and evapotranspiration (ET) at early (a), middle (b) and late (c) 
growing stages of maize 


5o: HERMAN; v: ERIN; no 冠 层 叶 面 边界 层 空 气动 力 阻力 ，r;: 地 表面 到 冠 层 面 的 空气 动力 阻力 ; nc 冠 层 面 高 度 到 参 


pm ~ = S C * c z . s 
考 面 高 度 间 的 空气 动力 阻力 。 T. : surface soil resistance; 7; : canopy resistance; 7, : aerodynamic resistance of boundary layer of leaf; 7, : 


" " a = " * 
aerodynamic resistance from surface to canopy; 7, : aerodynamic resistance from canopy to reference surface height. 


3 结论 与 讨论 与 此 同时 ， 本 文中 考虑 了 大 气 CO. 浓度 对 玉米 
SARA RABAT ME RMB, xl S cin dide Kel eg A Sl 
Rosas vue omeezanee, "l RECRAIIRTRBESAERMRUDE 
进而 估算 了 各 阶段 玉米 蒸腾 和 土壤 蒸发 并 利用 涡 。 实验 数据 对 CO 浓度 对 冠 层 阻 力 的 影响 做 进一步 修 
度 相 关 实测 数据 对 模型 进行 验证 ， 得 到 以 下 结论 E RA AMARAS FRAN, 仅 考虑 了 CO: 
1) 将 考虑 了 大 气 CO; 浓度 对 玉米 冠 层 影响 的 冠 浓度 变化 对 冠 层 阻力 的 影响 ,缺少 CO; 


层 阻 力 模型 耦合 到 双 源 的 Sw 模型 中 ,能够 更 为 精 。 “ 叶 面 积 指数 的 影响 ,这 将 在 后 续 工作 中 进行 完善 。 
准 地 模拟 玉米 整个 生育 期 不 同 生长 阶段 半 小 时 时 间 。 “参考 文献 References 
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